On the frequency band of the /-mode CFS instability 
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Rapidly rotating neutron stars can be unstable to the gravitational-wave-driven CFS mechanism if 
they have a neutral point in the spectrum of nonaxisymmetric /-modes. We investigate the frequen- 
cies of these modes in two sequences of uniformly rotating polytropes using nonlinear simulations 
in full general relativity, determine the approximate locations of the neutral points, and derive limits 
on the observable frequency band available to the instability in these sequences. We find that general 
relativity enhances the detectability of a CFS-unstable neutron star substantially, both by widening 
the instability window and enlarging the band into the optimal range for interferometric detectors 
like LIGO, VIRGO, and GEO-600. 
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I. INTRODUCTION 

The possibility of detecting gravitational waves from 
neutron stars is an exciting prospect for present and 
planned observatories, since it may be our only direct 
probe into the very core of the star (neutrinos being an- 
other possibility). In addition, neutron stars have a very 
unique spectrum of stable and unstable modes of oscil- 
lation which could be identified in gravitational wave 
observations, and indirect conclusions concerning the 
equation of state and internal structure of the star could 
be drawn similar to helio- and asteroseismology. 

Instabilities are particularly interesting in this context, 
since in this case a spectral line can be amplified either 
directly (by dynamical transition) or driven by radiative 
fields to potentially large saturation amplitudes, thereby 
enhancing the chances of detection. A number of neu- 
tron star instabilities are known (for a review see HI) 
which will be active subject to particular (different) re- 
quirements on the stellar characteristics. The resulting 
gravitational wave signals are expected to carry signa- 
tures unique to the unstable mode and the structure of 
the star. 

Typical pressure oscillations in neutron stars have fre- 
quencies in the kHz band, and are therefore, even if 
unstable, not located in the optimal range of sensitiv- 
ity of interferometric detectors like LIGO, VIRGO and 
GEO-600. However, a particular class of modes, the 
nonaxisymmetric /-modes, also have a counter-rotating 
branch which extends to very low frequencies in the 
range of hundreds of Hertz, which is optimal for de- 
tection by laser interferometers. These modes are the 
subject of the current study. 

These counter-rotating /-modes could grow to 
large amplitudes (and produce correspondingly larger 
gravitational-wave strains) in a manner first discovered 
by Chandrasekhar, Friedman and Schutz [2 ,3], in which 
the dynamical gravitational field itself is the agent driv- 



ing the instability. This so-called CFS mechanism oper- 
ates when a mode that is counter-rotating with respect 
to the star appears as co-rotating in the inertial frame, 
its frequency having crossed the neutral point (zero fre- 
quency in the inertial frame) due to the high rotation 
rate of the star. If the neutron star spins faster than crit- 
ical the instability could be active, with some qualifica- 
tions we will mention below. 

The most important one of these unstable modes is 
likely the lowest-order one with spherical harmonic in- 
dices I = m = 2. In Newtonian gravity, this mode be- 
comes unstable at spin rates which are too high for typ- 
ical models of (uniformly rotating) neutron stars. Gen- 
eral relativity enhances the instability substantially, as 
demonstrated first by Stergioulas and Friedman [4], and 
therefore makes it a more promising source for detec- 
tion. 

We will focus here on /-modes, which have pressure 
gradients as the restoring force. Since no successful lin- 
ear method exists to address this problem in a fully rel- 
ativistic stellar model, we employ full nonlinear numer- 
ical simulations to extract the mode frequencies. 

The CFS mechanism is not restricted to this class of 
modes. The only requirement is that, along a sequence 
of stars with increasing rotational over binding energy 
T / 1 W | , a mode that is counter-rotating in the rest frame 
of the star changes type to corotation in the frame of a 
faraway observer. Gravitational waves extract positive 
angular momentum from the star, and the mode's an- 
gular momentum becomes increasingly more negative 
with respect to the star, i.e. its amplitude increases with 
time, r-modes (a class of inertial modes) and ztf-modes 
(trapped spacetime oscillations) have been discussed in 
this context as well [5 ,0. 

The actual instability window depends on a number 
of parameters. The lower limit in terms of rotation is 
given by the neutral point, which for the £ = m — 2 
bar-mode in Newtonian gravity has been found to be 
close to T/|W| = 0.14. Higher order modes with I = m 
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have neutral points at lower T/|W|, i.e. they have a 
larger instability window. But the growth timescale of 
the instability also depends on m, such that higher or- 
der modes need a longer time to grow to high ampli- 
tudes. In that, the instability is competing with viscosity 
in the star which damps the mode amplitude. Instabil- 
ity requires 1/t = 1/t gw + l/T viscosity to be positive, 
where T denotes the growth timescale of the particular 
mechanism. In particular, shear and bulk viscosity in 
the neutron star provide damping agents which leave an 
effective temperature window around 10 8 . . . lO^K and 
mode numbers m < 5 171 . 

A superfluid core in the neutron star may also be in- 
strumental in suppressing the /-mode CFS instability 
by mutual friction between quantized neutron vortices 
||8j|9l. This can only occur below the transition tempera- 
ture to a superfluid state (s=a 10 9 i<C) and is therefore only 
relevant for older neutron stars. The effectiveness of 
mutual friction is based on estimates for the coupling 
strength. Also, current predictions for the instability 
windows assume Newtonian physics, while general rel- 
ativity is known to enhance the gravitational instability. 
Therefore it is not yet clear whether the instability could 
still be active also in some old neutron stars (the r-mode 
instability could also provide the dominant spin-limit 
mechanism in those cases). 

Another unknown is the saturation amplitude of the 
deformation. In the case of the CFS instability for the 
r-modes, previous studies have indicated that nonlinear 
mode-mode couplings appear to be effective in saturat- 
ing the instability at low amplitudes [10-13 1. A sim- 
ilar effect could occur for the /-mode instability, and 
whether this is the case or not is presently unknown. 
Since the spectrum of pressure modes is less dense than 
the inertia! mode spectrum and has a different charac- 
teristic eigenfunction, nonlinear mode-mode couplings 
could be less effective in this case. Also, some numeri- 
cal evidence for the development of the CFS instability 
in rapidly rotating Newtonian polytropes [14, 15 1 seems 
to indicate that large amplitudes could be reached (but 
these studies also enhance the growth rate by a large fac- 
tor, so they may not represent the actual mode coupling 
accurately). If that is the case, estimates using equilib- 
rium ellipsoids [16] predict that the gravitational wave 
signal would sweep through a band in the range of 100's 
of Hz and could be visible up to 140 Mpc in Advanced 
LIGO. 

The /-mode CFS instability requires a T/|W| of sev- 
eral percent, which is low compared to the dynamical 
bar-mode (T/ 1 W\ ~ 0.24) but still rather high compared 
to typical estimated iron core rotation rates. On the 
other hand, a subset of massive stars may undergo an 
evolution in which angular momentum transport from 
the core is less efficient than in typical cases, which 
is consistent with (but not supported by) the collap- 
sar model for long gamma-ray bursts which requires a 
dense disk of sufficient angular momentum as a source 
for powering the engine. To have a reasonable chance of 



detecting a CFS-unstable star, it is therefore important 
to reach high signal-to-noise ratios for a sufficient event 
rate. 

This study is focused on the frequency band and the 
location of the neutral point for the m = 2 and m = 3 
modes. The paper is organized as follows: in Section [fl] 
we describe the physical model, our numerical teciv 
niques, and the equilibrium stars used as initial data for 
time evolution. In Section III we present results from 
code tests and compare to previous work. Section |IV] 
contains the results of our study, and in Section [V] we 
conclude with a discussion. 



II. PHYSICAL MODEL AND NUMERICAL 
TECHNIQUES 

A. Basic physical model 

As mentioned in the introduction, we model neutron 
stars as (uniformly rotating) polytropes, i.e. those equi- 
librium stars whose structure follows from the pressure- 
density relation P = Kp T , with P the pressure, p the 
rest-mass density, and K and T free parameters. This 
model is simplified in a number of ways. Actual neu- 
tron stars are expected to consist of a composition of 
baryons and leptons which depends on the local ther- 
modynamical properties and the evolutionary history of 
the star IfTTI . In addition to protons and neutrons, the 
core of the neutron star may also contain hyperons or 
kaons, and it is possible that a deconfined equilibrium 
quark phase exists at high densities. These uncertain- 
ties in the microstructure affect the macroscopic struc- 
ture and dynamics as well, in particular since different 
compositions have a different stiffness in reaction to per- 
turbations. However, polytropes have the advantage to 
essentially only depend on one free parameter (the con- 
stant K can be viewed as a scaling parameter, since dif- 
ferent choices of K will lead to the same dynamics in a 
properly scaled system of units, see also the discussion 
in |4J), and since they have been used in many other 
studies they make it easier to isolate effects not associ- 
ated with the microstructure, e.g. the choice of Newto- 
nian vs. relativistic models. 

For the time development of perturbations we also 
assume the star to be an ideal fluid governed by the 
gamma-law P = (T — 1) p e, where we select T to be 
the same as in the polytropic relation, and e is the spe- 
cific internal energy density of a fluid element as mea- 
sured in a local rest frame. In addition, for numerical 
reasons, we restrict attention to purely adiabatic pertur- 
bations, i.e. we neglect local heat exchange by projecting 
the pressure back to P = Kp T - this procedure is com- 
monly referred to as assuming a polytropic equation of 
state. It is generally used to reduce artefacts associated 
with evolving the surface of the polytrope. 

In particular, these assumptions neglect non-ideal ef- 
fects like superfluidity [18] and solid crusts Ifl9] , which 



are important in neutron stars below 10 9 . . . 10 10 K, i.e. 
likely already a short time after birth. We also neglect 
magnetic fields by assumption of a pure ideal fluid. The 
resulting model is therefore not suited to give quantita- 
tively correct results, but allows us to qualitatively dis- 
cuss the effect of including full general relativity as part 
of the model, specifically in comparison to earlier results 
in the Cowling approximation ||201|2T1 . 



B. Numerical method 

We study neutron star oscillations by the time devel- 
opment of perturbed equilibrium stars in full general 
relativistic hydrodynamics. The numerical code used 
has been developed inside the Cactus computational in- 
frastructure [22J, and uses the Carpet mesh refinement 
and multi-mesh driver [23]. The module for evolving 
the general relativistic fluid fields is the Typhon multi- 
block code l24ll . which has been coupled to a multi- 
block implementation of the generalized harmonic form 
of Einstein's field equations [25 26 j. We point the reader 
to these publications concerning the details of the differ- 
ential equations and the discrete implementation. 

The generalized harmonic evolution code provides 
a linearly stable initial value problem by virtue of be- 
ing a first-order symmetric hyperbolic system. In ad- 
dition, we use maximally dissipative outer boundary 
conditions discretized with linearly stable penalty op- 
erators, and preserve the linear stability in a discrete 
sense by using one-sided finite difference operators sat- 
isfying the summation-by-parts property (see E5l l26l 
and references). The scheme also employs a constraint- 
damping technique and Kreiss-Oliger type dissipation 
operators of high order (we usually employ 8th-order 
accurate finite difference operators in the interior) to 
stabilize the non-linear system. These properties are 
particularly important for long-term stable evolutions 
of black holes and neutron stars. In particular, we are 
able to evolve isolated neutron stars for apparently arbi- 
trary long times with domain boundaries very close to 
the stellar surface, and see only minimal growth of con- 
straint violations. The accuracy of the numerical evolu- 
tion system is nowadays clearly dominated by the (low- 
order accurate) methods for evolving the fluid fields. 

For the fluid fields, we employ a standard finite- 
volume discretization with PPM reconstruction and 
HLL flux estimates. The evolution systems are cou- 
pled via the method of lines, and the time evolution 
is performed with 3rd-order Runge-Kutta integrators. 
Since the fluid scheme is unable to deal with matter- 
vacuum interfaces, we employ a low-density artificial 
atmosphere outside of the isolated star to regularize the 
method. 

Initial data for the time evolution is obtained with 
the rns code [27] and mapped to the numerical domain. 
Directly afterwards, a perturbation of the fluid fields 
is added to the equilibrium star to preferentially excite 



Model p c r p /r e M T/\W\ fi/n K 
xlCT 3 xlCT 2 
BUO 1.28 1.00 1.400 0.000 0.000 
BU7 1.28 0.65 1.666 8.439 0.850 



Table I: Properties of the models BUO and BU7 used for com- 
parison with previous publications. The columns are the cen- 
tral density p c , ratio of polar over equatorial coordinate radii 
r p I r e , gravitational mass M, ratio of kinetic over gravitational 
binding energy T/\W\, and ratio of angular velocity over the 
Kepler limit Cl/Ci K . All models are uniformly rotating poly- 
tropes with K = 100 and 1 = 2. 



the relevant oscillation modes. Mode frequency extrac- 
tion is performed by either extracting time profiles of 
fluid quantities at certain locations inside the star (in 
particular the center), or performing discrete Fourier 
transforms on coordinate circles in the equatorial plane 
Il28l l29t . The resulting time evolutions are then again 
Fourier-transformed in post-processing to obtain power 
spectra. 



C. Initial models 

As mentioned earlier, all initial models used in this 
paper have been generated with the rns code. A model 
is specified by its central density p c , the ratio of polar 
to equatorial coordinate radii, and the polytropic con- 
stants K and T. In addition, the resolution of the two- 
dimensional grid n s , n m can be adjusted, as well as the 
target accuracy for the iterative procedure e acc , and the 
maximal spectral coefficient £ max - We fix these values to 
n s = 601, n m = 301, e acc = 10~ 8 and i max = 10, which 
are sufficient to reproduce properties of published mod- 
els up to about 10~ 3 . . . 10~ 4 relative accuracy. 

For purposes of testing the code, we use models from 
the well-investigated BU sequence |20} I30H32I , which 
is a set of uniformly rotating polytropes with r = 2, 
K = 100 and p c = 1.28 x 10 -3 . 1 The properties of the 
nonrotating model BUO and the rapidly rotating BU7 
used for our comparisons are listed in Table II] For con- 



A note on units: The numerical code implicitly employs ge- 
ometrized units G = c = 1. Since polytropes are scale invariant, 
each choice of K needs to be supplemented by an additional choice 
of length to convert to cgs units. In fact, K = 1 would be sufficient, 
and we use different values only for numerical convenience. Any 
sequence of polytropes gives results for all possible models of the 
same stratification, i.e. the range of masses in terms of cgs units can 
be fixed by observational constraints. (A more precise analysis would 
make little sense here, given that polytropes are a substantial sim- 
plification of the actual interior structure of neutron stars.) For the 
sequence BU test models, we will fix the length scale to Mq = 1, 
which implies p c = 1.28 x Kr 3 c 6 /G 3 M| rs 7.906 x 10 14 g/cm 3 at a 
mass range M = 1.4 . . . 1.7M Q . 
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Model 



pc r v /r e 
xlCT 3 



M 



T/\W\ fl/Cl K 
xlQ- 2 



50 2.6812 1.0000 1.6287 0.000 0.000 

51 2.6812 0.9160 1.6801 2.001 0.380 

52 2.6812 0.8338 1.7367 4.001 0.552 

53 2.6812 0.7487 1.7995 6.001 0.704 

54 2.6812 0.6393 1.8700 8.000 0.894 



neutron star oscillations. For the case of stationary stars 
in the Cowling approximation without artificial pertur- 
bation, but covered by multiple blocks, results have al- 
ready been presented in [ 33 1 . In this study, we will not 
need the additional complexity of multi-block systems 
and instead resort to uniform grids to save computa- 
tional cost. 



Table II: Properties of the models SO to S4 used for comparison 
of the neutral point with Stergioulas & Friedman |4|. All mod- 
els are uniformly rotating polytropes with K = 100 and 1 = 2. 
The columns are the same as in Table [I] The particular mem- 
bers have been selected to give uniform spacing in T / | W| . 



Model 


Pc 


r p /r e M 


T/\W\ 


n/n K 




xl0~ 3 




xl0~ 2 




CO 


5.0 


1.0000 1.3056 


0.000 


0.000 


CI 


5.0 


0.8699 1.3598 


4.000 


0.472 


C2 


5.0 


0.7467 1.4499 


8.000 


0.686 


C3 


5.0 


0.6040 1.5609 


12.00 


0.910 



Table III: Properties of the models CO to C3 used for inves- 
tigating the effects of stiffness on nonaxisymmetric /-modes. 
All models are uniformly rotating polytropes with K = 1000 
and T = 2.5. The columns are the same as in TableU 



verting to cgs units later, we will assume Mq = 1. 

In addition to these models, we have constructed a se- 
quence of T = 2 polytropes with a higher central density 
p c = 2.6812 x 10 -3 , for purposes of comparing with the 
approximate estimates of the location of /-mode neutral 
points in full general relativity by Stergioulas & Fried- 
man |4|. This choice results in a central energy density 
of e c = 3.4 x 10~ 3 . The properties of selected models 
of this "S" sequence are listed in Table [TT] Finally, to in- 
vestigate the effects of stiffness in the equation of state, 
we also define a sequence with T = 2.5, which we will 
denote "C" sequence. Selected members are listed in Ta- 
ble III Note that the fastest rotating member, C3, could 
onlybe constructed by lowering the convergence factor 
of rns to 0.3. Since the mass of the nonrotating member 
is M(0 = 0) = 1.3056, we convert to cgs units with 
a length factor L C gs = 1.9226 x 10 5 cm: this results in a 
mass of the nonrotating neutron star of M(Q = 0) 
1.7 . V'l . , which is more consistent with expectations from 
observations. We repeat that this choice is arbitrary, and 
our results for frequencies can be rescaled to any other 
mass range as desired. 



A. Oscillations of a nonrotating neutron star 

Oscillation frequencies in nonrotating general rela- 
tivistic polytropes, and in particular in the BU0 model 
which we consider here, have been reported in l20l l3ll 
34 1 in the Cowling approximation (i.e. keeping the 
spacetime fixed) and in |32] in general relativity sub- 
ject to the conformal flatness condtion. To compare with 
these results, we have generated the model with rns 
and mapped it to a uniform grid with coordinate extent 
[— 13, 13] 3 . The stellar surface is located at coordinate 
radius R = 8.2, which means that the outer boundary 
is at ~ 1.6 times the stellar radius. This ratio is suffi- 
cient to avoid artefacts from surface-outer boundary in- 
teractions in the fluid system. The techniques used for 
the spacetime evolution, in particular the use of maxi- 
mally dissipative outer boundary conditions via penalty 
operators 2 and the constraint-damping property of the 
generalized harmonic evolution enforce stability of the 
spacetime system and only a moderate growth of con- 
straint violations. Experiments with larger coordinate 
domains (while keeping the resolution fixed) produce 
no difference in the resulting fluid mode frequencies 
within the given error bars. 

The spherically symmetric fundamental mode (usu- 
ally denoted F) is already excited by discretization er- 
rors, but to be consistent with the remaining study and 
with published work we impose a simple perturbation 
of the form (311 



Sp = Ap c cos 



nr 
2r~s 



(1) 



where Sp is the perturbation in rest-mass density, A an 
arbitrary amplitude, r the radial coordinate, and r$ is the 
coordinate radius of the surface. For this mode it is suffi- 
cient to extract the time evolution of the central density 
and perform a Fourier transform on the profile. Fig. n] 
compares the spectra obtained from the Cowling ap- 
proximation with fully general relativistic models. All 
data were obtained using a small perturbation ampli- 
tude (A = 10~ 3 ) and from an evolution over 10 ms. The 



III. CODE TESTS 

In this section, we present results from a number of 
standard code tests and compare to published results for 



2 In principle, constraint-preserving outer boundary conditions 1 35 
l36l could be a superior choice, but we find that for studies of 
fluid modes maximally dissipative boundary conditions give good 
enough results. 
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Figure 1: Fundamental mode in the spherically symmetric 
polytrope BU0 (r = 2.0, K - 100, central density p c = 
1.28 x 10~ 3 ): The mode is excited with a nodeless spherically 
symmetric density perturbation. This plot shows a discrete 
Fourier transform of the central density time evolution pc(t), 
for both a model in the Cowling approximation and in full 
general relativity. The resulting major peak is the fundamental 
mode F, which can be compared to published estimates. Since 
we do not impose an exact eigenfunction, the initial perturba- 
tion (and also the discretization error) also have some overlap 
with the first overtone HI, although the perturbation is node- 
less. The differenced plots refer to central finite differences of 
the Fourier transforms, which allow to locate the maxima of 
the peaks more accurately. (The difference graph values have 
been rescaled for illustration.) 



maximum evolution time we use is not limited by nu- 
merical stability, but by numerical dissipation present 
in the hydrodynamical scheme (this aspect will be dis- 
cussed in more detail in the context of the nonaxisym- 
metric mode extraction below). 

Fig. [T] shows that the fundamental mode frequencies 
are very different for models in the Cowling approxima- 
tion and fully general relativistic models (see also e.g. 
11321 ). In particular, we obtain a frequency of v(F) = 
2694(40) Hz in the Cowling approximation, which is in 
agreement with [31 J (2706 Hz), and for full general rela- 
tivity v(F) = 1465(40) Hz, compared to 1458 Hz in |32l . 
All error bars are rough estimates based on the distance 
between the root location in the first derivative of the 
Fourier spectrum and the next non-zero value. 

Non-radial modes in this model have been investi- 
gated as well, in particular the set of £ = 2 /-modes 
(which are degenerate in m since the background is 
spherically symmetric). To excite these modes, we im- 
pose a perturbation of the form 



5p = Ap c sin 



nr 
rs 



(2) 



and extract the value of the density at several locations 
between the center and the equator every few itera- 



tions. In post-processing, these time sequences are again 
Fourier-transformed to obtain the mode frequencies. In 
the Cowling approximation, we obtain, for the / = 2 
/-mode, a value of v( 2 f) = 1884(40) Hz, compared to 
1890 Hz reported in 1201 . In full general relativity we 
get v( 2 f) = 1601 (30) Hz, which is in marginal agree- 
ment with the result 1586 Hz from the CFC approxima- 
tion 153 . 



B. Oscillations of a rapidly rotating neutron star 

As a test case for a rapidly rotating model, we select 
the BU7 polytrope, which is defined by p c = 1.28 x 10 -3 , 
K — 100, r = 2 and a coordinate axis ratio r p /r e = 0.65. 
For the fundamental quasi-radial F mode, we select the 
same form of the perturbation function as above, but 
since the star is strongly rotationally deformed we mod- 
ify the surface coordinate rg to interpolate between the 
polar radius and the equatorial radius with an ellip- 
soidal approximation. Frequencies for this mode have 
been reported as v(F) = 2546 Hz in the Cowling approx- 
imation 11201 and v(F) = 1207Hz in the CFC approxi- 
mation to general relativity |32|. We obtain 2491(50) Hz 
and 1204(40) Hz in agreement with these values. 

For the axisymmetric £ = 2 /-mode, we again im- 
pose the same perturbation as in the nonrotating case, 
but with a ^-dependent value of rg, and extract mode 
frequencies of v( 2 f) = 1708(20) Hz in the Cowling ap- 
proximation, and v( 2 f) = 1694(35) Hz in full general 
relativity. The first value matches with 1703 Hz cited 
from [20 1, whereas the second agrees with the CFC ap- 
proximation (1720 Hz) 1521. 

Since we will be concerned with nonaxisymmetric /- 
modes later, we also investigate the £ = 2,\m\ = 2 f- 
modes in this model and compare to the recent study by 
Gaertig & Kokkotas [20J in the Cowling approximation. 
The techniques used here are also used to determine the 
nonaxisymmetric modes in full general relativity, and 
therefore we will describe the method in more detail. 

To extract fluid perturbations of a particular m, we use 
a method similar to the one in 11291 137|| : On the equa- 
torial plane of the star, we specify coordinate radii r ( , 
discretize the corresponding circles with sample points 
(usually 100 are sufficient), interpolate the density onto 
the sample points, and project the resulting density an- 
gular profiles by discrete Fourier transforms. The result- 
ing data output every few iterations gives a time profile 
of the Fourier amplitude of the nonaxisymmetric part of 
the density for each mode \m\ ^ in the analysis. We 
normalize these perturbations to the average density on 
the circle and denote the amplitudes A m to distinguish 
between the mode strengths. 

This method is particularly useful to identify the I = 
\m\ /-modes for the following reasons. First, the modes 
have a dominant pressure perturbation, which, by the 
adiabatic condition, also implies a dominant density 
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Figure 2: Example of a signal obtained by the mode extrac- 
tion technique described in the text. The background model is 
a member of the sequence C (number CI, see Table [Tn}, with 
K = 1000, T = 2.5 and p c = 5 x 10~ 3 . A nonaxisymmetric 
£ = 2, m = 2 density perturbation excites the corresponding 
two /-modes, and the resulting oscillations are extracted by a 
discrete Fourier analysis of the density on discrete rings. The 
plot shows the resulting signal of the mode Ai (for m = 2) for 
an evolution in full general relativity over 13 ms. The signal 
exhibits oscillations which are dominated by two frequencies. 
Also visible is a decay in the signal strength, which is a nu- 
merical artifact caused by the finite-volume scheme and the 
inaccurate discretization near the stellar surface. 

perturbation. Second, the spherical harmonics Ynji have 
meridional node lines <p = const, which give orthogonal 
projections with the corresponding Fourier basis func- 
tions in the equatorial plane. Third, the artificial stellar 
drift associated with the discretization error (which ap- 
pears as low-order components in the post-processing 
time Fourier transform of the axisymmetric average for 
m — 0) does not appear in the profiles for \m\ > 0. 

An additional tool we can use is the choice of a free 
phase function in the extraction. The mode projection 
on each circle is naturally written as a projection onto 
exp(im((p + <po)), where the phase angle </>o can be cho- 
sen freely. In particular, we will use a time-dependent 
function of the form <po(t) = Cit, where O is the ro- 
tational angular frequency of the stellar background 
model. This choice of phase shifts frequencies by a fac- 
tor of mCl, which is useful for extraction of mode fre- 
quencies near corotation (see also |20|). 

An example of a signal 3 obtained in this way is given 
in Fig. [2] The particular model used here is number 
CI of the C sequence (see Table [III), with K = 1000, 
T = 2.5 and p c = 5 x 10 -3 . The model is evolved, in the 



Cowling approximation, for about 13 ms. After an initial 
transient, the two frequencies corresponding to the two 
\m\ =2 /-modes are already visible in the signal. The 
signal also exhibits a decay over time. Gravitational ra- 
diation will reduce the mode amplitude over time (both 
I = \m\ = 2 /-modes are stable in this model), but on 
much longer timescales than in the plot. The decay is 
therefore a result of numerical dissipation in the finite- 
volume scheme, possibly in part also due to inaccuracies 
of the evolution near the surface. In practice, this decay 
is not a concern as long as a sufficiently large number 
of cycles are available to extract the mode frequency 4 , 
i.e. it limits the total time interval St available for mode 
extraction, and therefore also the total frequency resolu- 
tion in post-processing. 

There are three direct ways to reduce this effect: (i) In- 
crease the resolution. Since numerical dissipation scales 
with the grid spacing, a lower spacing also reduces the 
mode damping. In addition, the noise level is slightly 
lower when going to higher resolutions, and it is in 
part also determined by the (fixed) floating point accu- 
racy. Over all, higher spatial resolutions therefore also 
imply higher frequency resolutions in post-processing. 
In practice, three-dimensional simulations in full gen- 
eral relativity are still resource-intense, and given the 
systematical uncertainties in the choice of background 
model and microphysics we feel that a moderate fre- 
quency resolution as used here is sufficient, (ii) Increase 
the initial perturbation amplitude. Assuming the same 
damping, this will allow for a longer extraction win- 
dows and increase the frequency resolution. However, 
high initial amplitudes increase the relevance of nonlin- 
ear effects, which we do not intend to study here. Ex- 
periments show that this is probably not important until 
perturbation amplitudes larger than A = 0.01 are being 
used, but for safety we use A = 0.001 exclusively here, 
(iii) Use a different numerical scheme. Since the numeri- 
cal dissipation is a result of the hydrodynamical scheme 
(it also appears in evolutions in the Cowling approxima- 
tion) a change of scheme has the potential to improve 
the result. We have experimented with a monotonized 
central (MC) reconstruction instead of PPM, and found 
that with MC the damping is less prevalent. However, 
MC reconstruction also introduces a stronger angular 
momentum drift in the star. 

A direct Fourier transform of the mode signal from 
Fig. [2] in post-processing leads to the spectrum shown 
in Fig. [3] Both the co- and the counter-rotating / = 2 
/-modes are excited by the perturbation. As explained 
earlier, the frequencies are obtained from the roots in 
the derivative of the power spectral density, which we 



3 We will frequently refer to the numerically generated time series as 
a signal to underline that we are performing numerical experiments. 



4 Actually a damped oscillator has a different peak frequency than 
its undamped counterpart; we have estimated this effect from the 
damping timescale, but have found that it is far below the error bars 
reported in this publication. 
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Figure 3: Fourier transform of the signal from Fig. [5] obtained 
in post-processing. The initial I = 2, m = 2 perturbation ex- 
cites both the co- and counter-rotating /-modes. To determine 
the actual mode frequency, we do not use the location of the 
peak, but the root of the derivates obtained from second order 
central differencing, as shown in the inset. 



approximate with second order accurate centered finite 
differences. 

Previous results for nonaxisymmetric modes in 
rapidly rotating stars have been obtained, in the Cowl- 
ing approximation and using linearized equations, by 
Gaertig & Kokkotas [20J and Kriiger, Gaertig & Kokko- 
tas [21 1 . To compare with these results we have evolved 
the rapidly rotating model BU7 with our nonlinear code 
in the Cowling approximation. We obtain, for the coro- 
tating t = 2 /-mode in the laboratory frame, a frequency 
of 2447(25) Hz, compared to 2454 Hz in |20l, and for the 
counter-rotating mode 183(30) Hz, compared to 210 Hz. 

Overall, our code tests show that we can reliably re- 
produce published /-mode frequencies in rapidly rotat- 
ing stars both in the Cowling approximation and in dy- 
namical spacetimes. 



IV. RESULTS 

We have investigated the I — \m\ =2 and I — \m\ =3 
/-mode frequencies in two sequences of different stiff- 
ness, both in the Cowling approximation and in full gen- 
eral relativity. In the following two sections, we describe 
the results of our simulations for each sequence sepa- 
rately, and conclude this part with a comparison of the 
spectra in the two sequences. 



A. /-mode spectrum of sequence S (T = 2.0) 

The first set of polytropes is the S sequence as de- 
scribed in Table [TTJ. which have been chosen to com- 



Model T/\W\ 


Vrot H 2 f2) v{ 2 f-2) 




xl0~ 2 


Hz Hz Hz 


SO 


0.000 


2162 (20) 2162 (20) 


SI 


2.000 


644 1130 (35) 2888 (35) 


S2 


4.000 


889 630 (35) 3056 (30) 


S3 


6.000 


1059 215 (40) 3134 (30) 


S4 


8.000 


1183 -240 (30) 3154 (30) 



Table IV: Nonaxisymmetric £ = \m\ =2 /-mode frequen- 
cies in general relativistic polytropes in sequence S (Y = 2.0, 
p c = 2.6812 x 10~ 3 , see also Tablein). The quantity v ro t is the 
rotational frequency of the star, vp/z) is the frequency of the 
counter-rotating /-mode, and v( 2 /-2) the corresponding fre- 
quency of the corotating mode. The numbers in brackets are 
approximate error estimates. 



Model T/\W\ 


Vrot v c ( 2 f 2 ) V C ( 2 f- 2 ) 


A( 2 /2) 


A( 2 /-2) 




xl0~ 2 


Hz Hz Hz 


Hz 


Hz 


SO 


0.000 


2445 (35) 2445 (35) 


13.1% 


13.1% 


SI 


2.000 


644 1392 (35) 3176 (35) 


23.2% 


10.0% 


S2 


4.000 


889 873 (40) 3351 (40) 


38.6% 


97% 


S3 


6.000 


1059 427 (40) 3409 (35) 


98.6% 


8.8% 


S4 


8.000 


1183 -101 (55) 3382 (35) 


-57.9% 


7.2% 



Table V: Comparison of the Cowling approximation with full 
general relativity for the t = \m\ = 2 /-modes in the se- 
quence S. The frequencies obtained from the Cowling approx- 
imation are shown in the columns vc( 2 /2) an d vc( 2 f-2) (for 
the counter- and corotating modes). The columns Ap/2) and 
A( 2 /_2) display the relative error of the Cowling approxima- 
tion, defined as A = vq / v GR — 1. 



pare with the approximate location of the neutral point 
in general relativity by Stergioulas & Friedman |4J. 
The major surprise in that publication was that general 
relativity substantially enhances the CFS instability in 
comparison to calculations assuming a stellar structure 
given by Newtonian physics. However, the techniques 
used in |4| used a truncated gauge to determine the lo- 
cation of the neutral point, which should be less accu- 
rate particularly in the low-order modes most relevant 
for the CFS instability as an observable source. 

The I = \m\ =2 /-mode frequencies in models SO to 
S4 have been obtained with the Fourier extraction meth- 
ods described in the last section, and the resulting fre- 

As 



IV 



quencies (in the inertial frame) are listed in Table 
expected, the modes are degenerate in the nonrotating 
model and then encounter rotational splitting. The coro- 
tating modes saturate at some limiting frequency, here 
ps 3.2 kHz, whereas the counter-rotating modes have 
lower frequencies with higher stellar rotation, since they 
are dragged in the direction of the stellar rotation. 

To assess the influence of grid resolution on the ex- 
tracted frequencies, we have performed higher resolu- 
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Figure 4: Resolution test for the rapidly rotating model S3 (in 
full GR). Amplitude extraction results from the Fourier trans- 
form in both a simulation with 90 grid cells per coordinate di- 
rection and 120 cells are compared. The finer resolved simu- 
lation has also been run for longer time (up to 20 ms) to re- 
duce peak width and the resulting size of the error bars in the 
Fourier transform. The graph shows the Fourier transform, its 
derivatives and the null line to visually identify the location 
of the roots. It is apparent that the discrete peak maximum 
locations disagree slightly, but the locations of the root agree 
very well, better than the error bar estimate obtained from the 
half-width of the peak. 

tion simulations in selected cases (120 instead of 90 grid 
cells in each direction) and confirmed that the results are 
unchanged within the error bars. Fig.[4]shows an exam- 
ple of the result for the rapidly rotating model S3. 

The Cowling approximation is particularly inaccurate 
when low-order /-modes are considered [38 j, and there- 
fore we have also performed the same set of simula- 
tions with the spacetime variables kept fixed. The result- 
ing relative changes in the observed /-mode frequen- 
cies are detailed in Table |V] Modes obtained from the 
Cowling approximation generally show an overestima- 
tion of the actual frequency. In the nonrotating model, 
this amounts to a difference of about 13% in accord with 
earlier studies [32 J. For the corotating modes, the rela- 
tive error in frequency reduces when moving towards 
the mass-shedding limit to about 7%. However, in the 
counter-rotating modes, the error grows to rather large 
values as the sequence approaches the neutral point 
(and also beyond that). The definition of the relative er- 
ror is clearly singular at the neutral point; therefore, one 
could view this as a mathematical artifact of the defini- 
tion (in particular since the frequencies in the corotating 
frame will have a smaller error more consistent with the 
corotating modes). However, observationally the differ- 
ence between 427 Hz (Cowling) and 215 Hz (GR, both 
model S3) is very relevant, regardless of how the error is 
defined. 

The results from both the full simulations and those 



Figure 5: Frequencies of I = \m\ = 2 /-modes in the mod- 
els of sequence S (r = 2.0). The models are characterized by 
their ratio of rotational kinetic over gravitational binding en- 
ergy T/\W\, where T/\W\ = 0.08 is very close to the mass- 
shedding limit. Both results from simulations in the Cowling 
approximation (blue dashed graph) and those from full sim- 
ulations (black full graph) are displayed for comparison. The 
approximate location of the neutral point is indicated for each 
graph. 



using the Cowling approximation are displayed graphi- 
cally in Fig. [5] The shift between the Cowling approx- 
imation ana full simulations is almost constant over 
the sequence, which causes the large relative error at 
low mode frequencies. When parameterizing the se- 
quence with T / 1 W I as done here, the frequencies of the 
counter-rotating /-modes are almost a linear function, 
which makes a linear interpolation to locate the neutral 
point more reliable (this is in contrast to a parameteri- 
zation over the rotational frequency v, cf. ||20l ). In the 
Cowling approximation, the neutral point is located at 
T/\ W\ n = 0.076 ± 0.002, whereas in the full simulations 
we find T/\W\„ = 0.0695 ± 0.0018, i.e. the Cowling ap- 
proximation tends to underestimate the lower limit of 
the /-mode CFS instability window by about 10% in this 
particular case. 

The actual location of the neutral point is of some in- 
terest since we can compare it with the results obtained 
from perturbation theory. In [4], the location of the neu- 
tral point is quantified, for this particular sequence, at 
T/\W\ n = 0.0649 (Table 3 in [4]). This is outside of the 
approximate confidence level we have reported above 
by 3%, and in disagreement with the average value by 
5.9%. The most likely cause for this disagreement is the 
use of the truncated gauge in [4J since it is expected to 
be less accurate at lower order. We will provide some 
evidence for this possibility below. In this sense, the 
full numerical simulation seems to indicate a slightly 
smaller instability window than previously expected on 
grounds of the neutral point location. 
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Model T/\W\ 


Vrot 


f( 3 / 3 ) H 3 f- 3 ) 




xlO~ 2 


Hz 


Hz Hz 


SO 


0.000 





2758 (30) 2758 (30) 


SI 


2.000 


644 


1053 (30) 4067 (40) 


S2 


4.000 


889 


229 (30) 4405 (35) 


S3 


6.000 


1059 


-457 (40) 4535 (60) 


S4 


8.000 


1183 


-1307 (55) 4572 (50) 



Table VI: Frequencies of £ = \m\ =3 /-modes in the sequence 
S. Notations as in table Hvl 



Model T/\W\ Vrot V C ( 3 / 3 ) v c ( 3 f- 3 ) A( 3 / 3 ) A( 3 /_ 3 ) 





xl0~ 2 


Hz 


Hz Hz 


Hz 


Hz 


SO 


0.000 





2971 (60) 2971 (60) 


7.7% 


7.7% 


SI 


2.000 


644 


1222 (30) 4228 (30) 


16.1% 


4.0% 


S2 


4.000 


889 


382 (25) 4574 (25) 


66.8% 


3.8% 


S3 


6.000 


1059 


-331 (30) 4713 (35) 


-27.6% 


3.9% 


S4 


8.000 


1183 


-1232 (40) 4544 (40) 


-5.7% 


-0.6% 




T/IWI 



Figure 6: Comparison between £ = \m\ = 2 and £ = \m\ = 3 
/-mode frequencies along sequence S. The neutral point of the 
m = 3 mode is shifted with respect to the m = 2 mode. 



Table VII: Comparison of the Cowling approximation with full 
general relativity for the £ = \m\ =3 /-modes in the sequence 
S. Notations as in Table Ivl 



To investigate the influence of the mode order on the 
extracted frequencies, we also consider the £ = \m\ =3 
/-modes in the sequence S. The resulting frequencies are 
listed in TablefVlJ and the comparison between Co wling 
results and full simulations can be found in Table IVIli 
The frequencies of the corotating £ = 3 modes are gen- 
erally higher, in accord with expectations from pertur- 
bation theory: we have 2.8 to 4.6 kHz as opposed to 
3.2 kHz in the £ = 2 case. The Cowling approximation 
is more accurate for these modes. 

Fig.[6]shows a comparison between the £ = 2 and £ — 
3 frequencies obtained from the full simulations. The 
higher order mode becomes unstable at lower 5 T/|W|, 
hereatT/|W|„ = 0.0467 ± 0.0010. 

We again compare this value with predictions from 
linear theory: the neutral point for the £ = m = 3 /- 
mode from [4], Table 3 is stated as T/\W\ n = 0.0455, 
which is only in much better agreement with our full 
nonlinear result (2.6% when comparing average values). 
This behavior is expected from the use of a truncated 
gauge in 1 4 1, and therefore can be cited as evidence iden- 
tifying this approximation as the cause of the disagree- 
ment. In addition, the errors for the choice of a trun- 
cated gauge are estimated in [4] for a model with some- 



what smaller central density than used here. Using the 
results of |4] for the lower-density model (full gauge: 
T/\W\ n = 0.0463, truncated gauge: T/\W\ n = 0.0459) 
and assuming that the relative error between the full 
and truncated gauge remains roughly the same at the 
somewhat larger central density used here, we find that 
the disagreement between our current results and those 
published in |4J for the 1 = m = 3 neutral points is less 
than 1.7%, which is within our current error bars. 



B. /-mode spectrum of sequence C (T = 2.5) 



5 We remind the reader that these statements apply to ideal fluid poly- 
tropes, and that viscosity and superfluidity determine the actual in- 
stability window in real neutron stars. However, the neutral point 
is a lower limit for this window. 



The sequence C (cf. Table III is designed to investi- 
gate the influence of stiffness in the equation of state on 
the nonaxisymmetric /-modes and on the location of the 
neutral point. All models have a central density which 
is close to the central density of the maximum mass non- 
rotating model, similar to the S sequence. The main dif- 
ference is that sequence C admits substantially higher 
rotation rates compared to sequence S: while the mass- 
shedding limit lies around T/ 1 W| f» 0.09 in sequence S, 
sequence C has a limit beyond T/|W| « 0.12. This is 
relevant not only for the range of frequencies expected 
from the CFS instability, but also for the growth time of 
an unstable mode. 

We report the frequencies for the £ = \m\ =2 /-modes 
in gener al relativistic simulations of sequence C in Ta- 
ble VIII (note the comment about the choice of units in 
|II C[ > . The general behavior of the frequencies is compa- 
rable to the softer equation of state, but the lower com- 
pressibility of the material shifts the corotating mode 
frequencies into higher bands. In addition, the larger 
maximum rotation rate introduces a much higher limit- 
ing frequency of the counter-rotating /-mode, for which 
we find a lower limit already in the kHz regime. The dif- 
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Model T/\W\ 
xW- 2 


Vrot 

Hz 


y( 2 fi) y( 2 f-2) 

Hz Hz 


Model T/\W\ 
xlQ- 2 


Vmt 
Hz 


H 3 h) y( 3 f- 3 ) 

Hz Hz 


CO 0.000 





2527 (70) 2527 (70) 


CO 


0.000 





3282(50) 3282(50 


CI 4.000 


1046 


698 (50) 3668 (50) 


CI 


4.000 


1046 


332 (40) 5374 (40) 


C2 8.000 


1422 


-170 (50) 3894 (60) 


C2 


8.000 


1422 


-1104 (40) 5828 (50) 


C3 12.00 


1654 


-1029 (50) 3995 (60) 


C3 


12.00 


1654 


5864 (40) 



Table VIII: Frequencies of £ = | m\ = 2 /-modes in the sequence 
C. Notations as in table lVIIII 



Table IX: Frequencies of I = \m \ = 3 /-modes in the sequence 
C. Notations as in Table I Villi We were unable to extract to 
counter-rotating /-mode frequency in the most rapidly rotat- 
ing model, therefore it is left open here. 
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Figure 7: Influence of the equation of state on the £ = \m\ =2 
/-mode frequencies. The results are from simulations of se- 
quences S and C (cf. Tables [TT| and III also Fig. [5j. While the 
corotating modes directly show the effect of the lower com- 
pressibility of the T = 2.5 models in a shift in the frequency 
band, the counter-rotating modes have comparable frequen- 
cies in the admissible range of stellar modes, and consequently 
a similar location of the neutral point in terms of T/ — W — . 



ferences between full simulations and the Cowling ap- 
proximation are not listed here, but they show a similar 
shift as in the S sequence, which again becomes particu- 
larly important at lower frequencies. However, even for 
the fastest rotator C3 we have investigated, the Cowl- 
ing result is -7 28(45 ) Hz compared to -1029(50) Hz re- 



ported in Table 

-29% 



VIII 



resulting in a relative difference of 



In Fig. f7\ we compare the mode frequencies between 
the S andC sequences to illustrate the influence of the 
equation of state on the /-modes. The less compress- 
ible models show larger frequencies in the corotating 
modes, and extend to higher negative frequencies in the 
counter-rotating modes as observed above. The location 
of the neutral point in terms of T/ 1 W| is very similar 
when comparing these particular sequences. 

The i = | m| = 3 /-modes in sequence C are listed 



in Table IX We were unable to extract the frequency of 
the counter-rotating mode in the most rapidly rotating 
model C3 due to numerical noise levels, therefore this 
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Figure 8: Same as Fig.[7j but for the I = \m\ =3 modes. 



information is left out of the table. We display a com- 
parison to the same modes in the sequence S in Fig. |§J 
The neutral point is again only mildly different between 
these sequences, as are the frequencies of the counter- 
rotating modes. 



V. CONCLUSION 

We have analyzed the spectrum of £ = \m\ =2 and 
I = \m\ =3 /-modes in rapidly rotating neutron stars in 
full general relativity. As mentioned in the introduction, 
these modes are interesting not only for general pur- 
poses of neutron star asteroseismology, but particularly 
because they have a lower frequency (counter-rotating) 
branch, and since a part of that branch can be unstable 
to the CFS instability. Our study allows us to investigate 
the observable frequency band covered by the CFS in- 
stability (if it is active) within the restrictions imposed 
by using a polytropic model. 

The location of the frequency bands available to the 
CFS /-mode instability is displayed in Fig. [9] The di- 
agram also shows the noise curves of Advanced LIGO 
and the proposed Einstein Telescope to compare with 
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Figure 9: Frequency band of the CFS instability based on the 
extraction of m = 2 and m = 3 /-mode frequencies in two 
sequences of uniformly rotating polytropes. Each band is con- 
strained by the neutral point (low frequencies on the left end) 
and the frequency of the counter-rotating mode in the most 
rapidly rotating model. For the m = 3 mode in the T = 2.5 se- 
quence of polytropes, the counter-rotating mode frequency is 
extrapolated from figure|8]to around « 2.5kHz (a lower limit 
is given by 1.1 kHz). To compare with the frequency varia- 
tion of the detector sensitivity, the noise curves (in 1 / \/Hz) of 
Advanced LIGO and the proposed Einstein Telescope are also 
displayed. 



the characteristic frequency dependence of the detector 
sensitivity. Actual strain amplitudes of the CFS instabil- 
ity can of course not be determined with the results of 
this study, but it is instructive to see how the variation 
of the mode order and the stiffness in the polytropic EOS 
affects the available frequencies. 

The band is determined by the location of the neu- 
tral point (where the pattern speed introduces very low 
frequency gravitational waves) and the frequency of the 
counter-rotating mode of the fastest rotator. If we take 
the approximate development of the instability as de- 
scribed by Lai and Shapiro (see figures 4 and 6 in 1 16 ]) at 
face value, we expect the signal of an unstable star to rise 
to high amplitudes with frequencies close to the nonax- 
isymmetric mode which caused the instability and then 
sweep through the band down near the neutral point to 
a stable state in which the star has lost its angular mo- 
mentum due to gravitational radiation. However, since 
the modes have to compete with viscosity, and since 
the growth time is shorter for stars with rotation rates 
which are farther removed (for that particular mode) 
from the neutral point, it appears more likely that mod- 
els near the breakup limit are relevant for observation 
than those near the neutral point. If that is true one 
would expect the full band displayed here to be a fair 
assessment of the actual observable frequency sweep in 
these particular models. It goes without saying that the 



actual band also depends on the particulars of different 
polytropic sequences (different central density), scales 
trivially with the polytrope mass range, and will be dif- 
ferent when more realistic stellar models are assumed. 
However, since the low-order /-modes mostly depend 
on average density and compressibility of the star, it is 
quite possible that the systematic scalings observed in 
these bands readily transfer to more realistic models. 

From Fig. [9] we can see that the upper third of the 
m — 2 frequency band of the T = 2.0 polytropes in- 
vestigated by Stergioulas and Friedman |4J lies in the 
optimal sensitivity range of both Advanced LIGO and 
the Einsteine Telescope. As explained above, this part of 
the band is also expected to be the most relevant one for 
purposes of detecting the CFS instability. We can also 
see from the diagram that the Cowling approximation 
is particularly inaccurate in this case, since it underesti- 
mates the upper band by over 60% (~ 100 Hz in model 
S4 as opposed to the full evolution value of ~ 240 Hz). 
This disagreement, however, diminishes in sequences 
with a larger ratio of the Kepler frequency to the rota- 
tional frequency of the neutral point. 

When increasing the mode order to m = 3, the band 
extends to about 1 .3 kHz, such that the most likely ob- 
servations (from the larger growth time of the instability 
farther from the neutral point as well as the larger num- 
ber of cycles, see |16|) will occur in a band around or 
below 1 kHz. A similar change occurs when changing 
the sequence to a much stiffer polytrope with T = 2.5, 
but retaining the mode order to quadrupole. While the 
emission of the most unstable polytropes is not any- 
more in the optimal regime of around 100 — 200 Hz it 
is still expected below 1 kHz. Finally, the full extent of 
the m — 3 band in the T = 2.5 sequence could not 
be deduced from our data due to numerical problems, 
we know a lower limit of around 1.1 kHz and, by ex- 
trapolation in the mostly linear part of the l/(T/|W|) 
diagram (cf. Fig. [§}, should be around 2.5 kHz. From 
earlier work on the location of the neutral point when 
changing the central density of the star [4J it is expected 
that sequences whose nonrotating members are further 
from the fundamental mode stability limit dM/ dp c = 
should have a neutral point at higher ratios of the crit- 
ical to maximal rotational frequency, and therefore will 
be more stable to the CFS mechanism, and would also 
occupy a smaller frequency band of instability. 

With the present study, we have investigated the ac- 
tual mode frequencies of particular nonaxisymmetric 
modes in uniformly rotating, general relativistic poly- 
tropes, and were able to infer information about the 
general location of the frequencies the instability could 
cover. There are a number of important directions for fu- 
ture work: The first one concerns the rotational profile of 
the star, since differential rotation opens up a larger in- 
stability window by the simple fact that it admits T/ \ W\ 
beyond the Kepler limit for uniformly rotating stars, and 
it is expected that proto-neutron stars are differentially 
rotating. For recent work on nonaxisymmetric modes 
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in rapidly differentially rotating stars in the Cowling 
approximation see |21J. The second limitation is the 
restriction to a polytropic stratification. A number of 
equations of state for cold neutron stars are available 
and can be used to more closely approximate the struc- 
ture of a young neutron star. Alternatively actual struc- 
tures of proto-neutron stars obtained from core-collapse 
simulations could serve as a guide to construct more re- 
alistic models. The third uncertainty concerns the role 
of superfluidity in modifying the mode spectrum (and 
thus both the location of the neutral point and the fre- 
quency band). In addition, mutual friction in a super- 
fluid core may modify the instability window in older 
neutron stars or even suppress the instability entirely. 
Related to this, it is important to gain a better under- 
standing of the role of bulk and shear viscosity in af- 
fecting the instability window and the frequency band 
of the /-mode CFS instability. Finally, it is important to 
understand the nonlinear couplings that will lead to a 
saturation of the mode amplitude as well as the over- 
all nonlinear development of instability, including the 
time-evolution of the expected gravitational-wave sig- 



nal. 
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